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Abstract 

We present a model to describe an underdamped long Josephson junction coupled to a single- 
mode electromagnetic cavity, and carry out numerical calculations using this model in various 
regimes. The coupling may occur through either the electric or the magnetic field of the cavity 
mode. When a current is injected into the junction, we find that the time-averaged voltage exhibits 
self- induced resonant steps due to coupling between the current in the junction and the electric field 
of the cavity mode. These steps are similar to those observed and calculated in small Josephson 
junctions. When a soliton is present in the junction (corresponding to a quantum of magnetic flux 
parallel to the junction plates), the SIRS's disappear if the electric field in the cavity is spatially 
uniform. If the cavity mode has a spatially varying electric field, there is a strong coupling between 
the soliton and the cavity mode. This coupling causes the soliton to become phase-locked to the 
cavity mode, and produces step-like anomalies on the soliton branch of the IV characteristics. If 
the coupling is strong enough, the frequency of the cavity mode is greatly red-shifted from its 
uncoupled value. We present simple geometrical arguments which account for this behavior. 
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I. INTRODUCTION 

Josephson junction arrays have now been studied for more than twenty yearsj^. They are 
of interest for many reasons. On a fundamental level, many unusual physical phenomena 
may occur in them (Kosterlitz-Thouless transitions, and quantum phase transitions from 
superconductor to insulator). On an applied level, they have been investigated as possible 
sources of coherent microwave radiation, which would occur if all the junctions in the array 
were to radiate in phase. 

These studies have also been extended to Josephson junctions, and arrays of junctions, 
coupled to a resonant cavity. When such a system is driven by an applied current, it may 
exhibit equally spaced self-induced resonant steps (SIRS's) at voltages equal to multiples 
of hfl/2e, where Q is the cavity resonance frequency |i4i]. More recently, it was shown that 
when an array of underdamped junctions is placed in a single-mode cavity, and the array is 
biased on one of the SIRS's, the array phase- locks and radiates energy coherently into the 
cavity, provided that the number of junctions exceeds a critical threshold^, 0, Q, 0. The 
SIRS's are thought to occur because the cavity acts back on the array like an a. c. driving 
field, producing voltage plateaus analogous to Shapiro steps. The physics of the threshold is 
believed to be that, since each junction is coupled to the same mode, it is effectively coupled 
to all the other junctions. Thus, there is a long range interaction for which phase locking 
occurs when the coupling strength, multiplied by the junction number, exceeds a critical 
value. 

Several models have been developed to describe these steps jo], [l^ These 
models generally succeed in reproducing many of the salient features of the experiments: 
SIRS's at the expected voltage, and an increase in the energy in the cavity proportional to the 
square of the number of junctions. The models have now been extended to two-dimensional 
(2D) junction arraysQ], where they show that, on a given step, the 2D a rray s radiate 



i rray g 



much more energy into the cavity than the ID arrays. However, recent work[l^ indicates 
that, in contrast to experiment, the numerical models always result in the junctions' being 
phase-locked to each other whenever they are biased on a SIRS. 



These theories have also been extended into the quantum regime |l5l. Il6l. Il7lll8lll9ll20l |21 1 . 
This regime is relevant when the junctions are very small, so that the non-commutation of 
the Cooper pair number operator and the phase operator is important. In this regime. 
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the states of the junction (or junctions) and the cavity become entangled, and uniquely 
quantum phenomena such as Rabi oscillations may be detectable in suitable experiments. 
Since such junction-cavity systems may be controllable externally, they may be useful as 
qubits in quantum computation. For this reason, the quantum regime has lately received 
extensive attention. 

Several groups have also investigated the IV characteristics of a Josephson junction cou- 
pled to a phonon mode. This system is formally analogous to the junction-resonant cavity 
system, in that both involve coupling between a junction and a harmonic oscillator mode. 
Helm et a/.j^, have developed a model describing the coupling between longitudinal 
phonons and intrinsic Josephson oscillations in cuprate superconductors, which, when suf- 
ficiently anisotropic, behave like stacks of underdamped Josephson junctions^, 25. l26|. 
These theories have also been extended into the quantum regime [27!. 

In all the above calculations, each superconducting island was treated as a sma// object, 
with only two degrees of freedom: the phase of the superconducting order parameter, and 
the number of Cooper pairs on each island. If the island is still large enough for Cooper pair 
number and phase to commute, then these two degrees of freedom can be treated classically. 
One then obtains a set of classical coupled second-order differential equations describing the 
dynamics of the junction array. For small islands, the phase and number are non-commuting 
operators, but there are still only two such operators per island. 

But one may often be interested in /on^ junctions, and arrays of long junctions. In such 
systems, the phase difference across the junction typically depends on position along the 
junction. When such a junction is coupled to a cavity, its behavior may therefore differ from 
that of small junctions. 

In this paper, we present a simple model for a long Josephson junction coupled to a single- 
mode resonant cavity. We present our model in the classical regime, in which number and 
phase degrees of freedom commute. However, the model should be readily extendible into 
the quantum regime. We consider both spatially uniform and spatially varying coupling 
between the junction and the cavity. If the coupling is uniform, we find that, as in the 
case of small junctions, there are SIRS's in the IV characteristics provided that there is no 
magnetic field parallel to the junction plates. But when such a magnetic field is introduced 
(in the form of a soliton propagating in the junction), the SIRS's are no longer present for 
uniform coupling. We will give simple arguments why the steps are not to be expected in 
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this case. If the couphng is non-uniform, however, we find that a moving sohton in the long 
junction does couple to the cavity mode. This coupling perturbs the cur rent- volt age (IV) 
characteristics of the junction, producing step-like structures in the IV characteristics, and 
excitation of the cavity mode. If the coupling is strong enough, this interaction produces a 
substantial shift in the frequency of the cavity mode. All these perturbations and steps can 
also be understood by simple physical arguments. 

As in the case of short junctions, some previous work has also been done treating the 
coupling of intrinsic long junctions to phonons. In particular, Preis et a/.j3] have treated 
this coupling in systems with a magnetic field parallel to the junction plates. They also 
find resonances and phonon frequency shifts analogous to those found in the present work, 
although they use a somewhat different mathematical model. Besides this work, several other 
groups have studied the interaction of solitons in long Josephson junctions with various types 



of harmonic modes. For example, Kleiner 



29| has investigated the interaction of solitons in 



stacks of long junctions with the cavity resonances whose frequency and wave number are 
determined by the stack geometry and by the Josephson plasma frequency. Machida et 
ffl/. jsfljl have numerically studied emission of electromagnetic waves by solitons in stacks of 



long i unctions which are coup, 
a/.bll, and Gronbech- Jensen 



ed to the Josephson plasma resonances of the stack. Salerno et 
, have considered the coupling of solitons in long junctions 
to externally applied microwave fields. Below we compare some of these results with our 
own. 

Since SIRS's may occur in long junctions, they may be useful as sources of sub-THz 
radiation. Hence, stacks of long junctions may be even more useful as such sources. Since 
the most anisotropic high-T^ cuprate superconductors behave as stacks of long underdamped 
Josephson junctions, these natural materials, too, may serve as a source of coherent sub-THz 
radiationj33|. 

The remainder of this paper is arranged as follows. In Section II, we derive the equations 
of motion for treating a long junction coupled to a single-mode resonant cavity. In Section 
HI, we present numerical results obtained from these equations of motion, for both no- 
soliton and soliton initial conditions, and for both uniform and non-uniform junction-cavity 
coupling. Finally, in Section IV, we discuss the numerical results and their experimental 
implications, and also suggest some possible extensions of this work. A detailed derivation 
of the Lagrangian for the junction-cavity interaction is presented in an Appendix. 
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II. DERIVATION OF EQUATIONS OF MOTION. 



A. Lagrangian for a Long Josephson Junction 

We will consider a long Josephson junction characterized by a gauge-invariant phase 
difference (f){x,y,t), where x and y are spatial coordinates and t is the time. (f){x,y,t) 
represents the phase difference between the two plates of the junction at points {x,y). We 
assume that the junction has length in the x direction and Ly ^ in the y direction. 
For mathematical convenience, we assume periodic boundary conditions in the x direction, 
so that 0(x + Lx,y,t) — (f){x,y,t). We also assume that the phase varies only in the x 
direction. 

In the absence of the cavity, the junction Lagrangian is Ljunc — I dxCjuncLy The La- 
grangian density jCjunc — K — U, where K and U are the kinetic and potential energies per 
unit area. K = edE'^ / {d>T^) is the energy density stored in the cavity electric field, d is the 
spacing between the junction plates, and e is the dielectric constant of the material within 
the junction. 

With the help of the Josephson relation 0^ — ^ — where V is the voltage drop 
across the junction, e is the magnitude of the electronic charge, and the subscript denotes a 
time derivative, this becomes K — g^^i^- 

The potential energy density U = Uj + Ub- Uj = cos(j) is the Josephson coupling 
energy density, Jc being the critical current density. The magnetic field energy per unit 
junction area is 

where B = V x A is the local magnetic induction. The vector potential A is related to the 
gauge-invariant phase by (f){x,y) = 6{x,y) — / A{x,y) ■ dl, where 6 is the local phase 
difference across the junction in a particular gauge, $o = hc/{2e) is the flux quantum, and 
the second term is the integral of the vector potential across the junction at the point (x, y). 
We now assume that B = By, i. e., is parallel to the plates of the junction. If we choose 
the gauge A = —A{x)z, then y) — 9{x, y) -\- ^A{x)d and therefore 

where the subscript denotes a partial derivative with respect to x. 
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If the superconductor in the upper plate has a complex order parameter 
exp{i6u), then the corresponding in-plane current density may be written 



2m* 



m*c 



-A|^, 



(3) 



where e* = 2e, m* is the effective mass of a Cooper pair, and the gradient is taken in the xy 
plane. The coefficient of A in eq. (jHl) can be identified with the London penetration depth 
A of the material in the upper plate, via the relation 47r(e*)^|\E'up/(m*c^) = Similar 
equations hold for 9i on the lower plate. Hence V^^ = V(0„ — 9i) satisfies 



ve 



(4) 



where we have assumed that materials in the upper and lower plates have the same pene- 
tration depth. Substituting eq. (H)) into eqs. (HJ and ((H), we obtain 



U. 



B 



d 

8^ 



0(Px 



{Ju — Ji, 



27rd dc 

where we have assumed that (f) varies only in the x direction, as is reasonable for B||y. 
Given Cjunc, the equations of motion can be obtained from 



(5) 



9 dCj^^Q 
dt d6f 



+ 



dx 



0. 



(6) 



In obtaining these equations and henceforth, we assume J„ — = 0. The equations then 
take the form 



1 



1 



sm ( 







(7) 



where (f = c^/e and 1/Aj = Aire* dJc/{hc^) is the squared Josephson penetration depth. 
This is the well-known equation for the phase in a long Josephson junction, in the limit 
of no dissipation. It leads to many solutions, including the sine-Gordon soliton JssJ 



0(x,t) 



4 tan 



exp 



{x-xo)/Xj~f3t 



where P = v/c is a scaled velocity. This corresponds 



to an excitation traveling in the positive x direction with speed v. This excitation carries 
a single quantum, $0 = hc/2e of magnetic flux in the positive z direction. The soliton 
is relativistically Lorentz-contracted and cannot exceed a speed of c. A similar antisoliton 
travels with constant speed v in the negative x direction. 



6 



B. Cavity Lagrangian 

We assume that the cavity supports one harmonic oscillator mode, described by a "dis- 
placement" variable and its time derivative g^- A suitable Lagrangian for this mode is 
Lose = ~ where M is the "mass" of the oscillator mode and K is the "spring 
constant." The corresponding Lagrange equation of motion is 



which gives gV + ^'^qr = 0, where Vt = (K/My^'^ is the oscillator frequency. Here, we 
envision the oscillator as an electromagnetic mode of a suitable resonant cavity. In this case, 
for some types of electromagnetic modes, qr is proportional to the electric field of the mode. 
The formalism described here would, however, also apply to suitable single-mode mechanical 
oscillators. With some modifications, it would also apply to electromagnetic modes in which 
qr represents the magnetic, rather than the electric field. 

C. Cavity-Oscillator Coupling and Bias Current 

We assume a capacitive coupling between the junction and the oscillator of the form as 
is shown in the Appendix, the coupling Lagrangian takes the form 



This is a natural extension of coupling assumed in Ref. jllj to a long Josephson junction 
[see, in particular, eq. (38) of that paper]. If the cavity electric field is non-uniform, the 
coupling depends on position along the cavity. In the numerical examples given below, we 
consider both position-independent and position-dependent coupling. 

Besides the capacitive coupling, there could, in principle, also be an inductive coupling 
between the junction and the cavity mode. In the present work, we do not include this term. 
However, if the cavity electric field has a non-zero curl, there is a corresponding magnetic 
field which is already included in eq. Q . We include this type of electric field in some of our 
calculations below. A more general derivation of the junction-cavity coupling can be found 
in the appendix. 

We also need to include a term in the Lagrangian corresponding to the bias current. This 
takes the form Lcurr = ^ / dx(p{x)Ly. 




(8) 
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D. Equations of Motion 



The equations of motion are obtained from the analogs of eqs. © and (jHI), but using 
the full Lagrangian. We derive the equations of motion assuming a sinusoidal coupling 
9e{x) = gEsm{k x), where qe is a constant. The equations of motion for uniform coupling 
are discussed below. 

The total Lagrangian is 



Ltot J i^^junc ^coup) dxLy -\- Lose -^curr- 



The Lagrange equations of motion take the form 



d_ 
di 



d_ 

dx 
d 



d 



i^^junc ~l~ ^ coup J 



(10) 



coup J 



(11) 



and 



d (dL 



tot 



dL 



tot 



0. 



dt \ dqr J dqr 
Carrying out these operations, we obtain the equations of motion as 



- (^sin0 - + g' sm{kx)qr = 0, 



where 



and 



2^2 



(12) 



(13) 



(14) 



(15) 



47rc/(e*)2M^ J —y-^^— v 
The above equations of motion are derived for a cavity electric field which varies sinu- 
soidally with position. To obtain the equations of motion for a uniform cavity electric field, 
one simply replaces the term g' sm.{kx) in eqs. (fT^ and (fT3|) (and subsequent equations) by 

g'. 

Eqs. (fT^ and (fT3j) do not include any damping. We simply incorporate damping by hand, 
by adding the appropriate terms to these equations. The resulting equations take the form 



1 



''XX —T) ' 



QjC 



n2 



J. 



Jr 



— ^ ( sin r ) + fi'' sin(A; a;)gV = 0, 



(16) 
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and 

Qr + g-gr + ^'^Qr = 4^^(g.)2 j^f ^ 7 sm{k x)(pttdxLy. (17) 

Here we have introduced dimensionless junction and cavity quality factors Qj and Qc, and 

a Josephson plasma frequency ujp = c/Xj = [ATre*dJc/ {he)]^^"^ . In terms of the junction 

parameters, Qj = eujjj cr, where a is the conductivity of the medium within the junction. The 

additional terms in eqs. (fT^ and (fTTj) insure that, in the absence of cavity-junction coupling, 

the equations of motion reduce to the standard results for a long Josephson junction with 
I I 

damping|35|], and for a damped harmonic oscillator. 



E. Reduction to a Set of Coupled First-Order Equations 

Eqs. p6|) and ()17|) are conveniently solved numerically if they are converted into a set 
of coupled first-order differential equations. Thus, we introduce the momenta canonically 
conjugate to and 0, namely Pr = Mq^ — Eqq' J sm{k x)(f)tdxLy and = Eocpt/c^ — 
g' sm{k x)Eoqr, where Eq = ji^^f^y- Substituting these variables into eqs. (fTH|) and (fTTj) . 
and rearranging, we obtain the following set of four coupled partial differential equations: 



K \M 



+ j sm{k x)p^dxLy 



Pr = -Mil Qr - -J7-qr 

(l>t = —P4, + g'sm{kx)c^qr 
Eq 

EqujI ( Jz 



{p^)t = Eo(j)xx ^ I sin ^ 

ujpEo (j)t 



(18) 



where K = l- / sin^{kx)dxLy 



F. Equations of Motion in Dimensionless Form 

It is convenient to rewrite these equations using a dimensionless time r = ujpt, length 
^ = x/\j, and cavity frequency Vl = Vt/up. We also introduce the dimensionless variables 
Pr = Pr/iMc), p^ = {Mc^/E^fip^/iMuOp). Eo = Eo/{Mc^), = g,/Aj, and finally 
g = g'LyC?. Finally, we require that k = 27im/Lx, where m is an integer, consistent with 
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periodic boundary conditions. With these substitutions, and further rearrangement, the 
equations of motion become 



EoP^ + g-rir sin( — - — , 

L/y 

EoVj, = (hff - [ sin (h '^''^ ^ 



JcJ Q 



J 



Pr = -^Qr-^^Qr, (19) 

where the dot is a derivative with respect to r and K = 1 — J sin^(/ca;) dx. 



III. NUMERICAL RESULTS 



A. Algorithm 

We solve eqs. ()19|) numerically by discretizing them on a spatial scale of A. Thus, 
the long junction. In discretized form, eqs. |T9|l become 



.j+i zfflj^ffl., 1 ^ ^ffiiQYQ j = ,,,^N, and = is the number of discrete sections of 



1 

'k 



N o 

+ A^o^Z^P^. sm( 



0i = Eop^^+g— sm(^^)g^ 
^ ~ (pi+i - + (f)i_i / 



Qj 

Pr = -^hr-^^Qr. (20) 

With this discretization, we are basically treating the long junction as N inductively coupled 
small junctions, each of which has critical current Ic = Jc{LxLy). 

We have solved eqs. (j2(Jj) numerically using a constant-time-step fourth-order Runge- 
Kutta method with a time step Ar = 0.001. We begin the simulation by initializing the 
variables pr^ Qr, p(t>j and the parameter ^ to zero; we have made various choices for initial 
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FIG. 1: IV curve for single long junction, with Qj = 10.0, A = 0.05, Q. = 2.5, -EgA^ = 1.0 x 10"^, 
Qc = 10 and all the phases initialized to zero. A SIRS is located at < y > /{NRIc) ~ vr, where Ic 
is the critical current of one of the individual small junctions. 
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FIG. 2: V{t)/{NRIc) versus dimensionless time r for small junctions no. (a) 59 and (b) 99 for 
4000 < T < 4010, plotted for a voltage on the SIRS of Fig. [H 



values of as discussed below. For a given Jz/Jc, we integrate the differential equations 
from r = 0tor = 5x 10^, then evaluate the voltages by averaging over the last r = 2 x 10^ 
units of time. The ratio 4^ is then increased or decreased by 0.01 and the set of equations 

•J c 

is solved again. In all cases, we use periodic boundary conditions in the x direction. 

We have considered two different initial phase distributions. To model a junction contain- 
ing no soliton, we choose all the initial phases 0j = for j = 1, A^. To model a junction 
containing a single soliton, we assume an initial configuration 0„(r = 0) = 2(| + tanh ^^), 
for n = 1, A^. As has been discussed in Ref. 2^1, for example, these boundary conditions 
are consistent with the presence of a single soliton, which carries one flux quantum. Our 
numerical results, as presented below, do show evidence for a soliton with these boundary 
conditions. 
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FIG. 3: Cavity variable qr{T) for the long junction in Fig. ^ at an applied current on the SIRS, 
plotted versus time r for the time interval 4000 < r < 4010. 



B. Results with No Soliton 

In Fig.0 we show the IV curve for a single long junction coupled to a cavity, using the no- 
soliton initial conditions (all phases chosen to equal zero) . We choose a position independent 
junction-cavity coupling with a constant qEqA = 1.0 x 10^^ and = 120. We plot the 
time- averaged voltages {V) in units of NRIc, where Ic = JdL^Ly) is the critical current 
of a single small junction. The IV characteristics have a step at {V)/{NRIc) = AttCi/Qj. 
This is a SIRS, similar to that seen in individual small junctions for a similar model (see 
Ref. (ll ), and occurs at {V) / {N RIc) = vr. This is the voltage expected from the model 
of Ref. Uj, which predicts that the SIRS's will occur when is an integer multiple of 
47r^. The step in Fig. ^occurs at 47?^^ = vr for n = 1. On this step, the phases of all 
the junctions oscillate coherently and are locked onto the cavity mode. To illustrate this 
coherence, we show in Figs. Efa) and|H^b) plots of the voltages for small junctions nos. 
59 and 99 over the same time interval. The plots show that the two junctions are indeed 
oscillating periodically and in phase with one another. In fact, the plots for all the 
small junctions are identical on this step. 

The phase locking between the junctions and the cavity is illustrated in Fig. El This 
Figure shows that the period of oscillation of qr{T) is identical to those of shown in 



Figs. 121 both equal ^ 



NRIc 

2tt 
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C. Results with Soliton Present. 

The results shown in Figs. are very similar to those obtained in Refs and js^ for 
a single small junction coupled to a cavity. However, we see very different results if there is 
a soliton initially present in the junction. In this case, it is important that the equation of 
motion for the long junction be discretized on a fine scale. If is too small, there are spurious 
steps in the IV characteristic produced by locking of the soliton to certain linear excitations 
generated purely by the numerical discretization. These steps, which are numerical artifacts 
of a too coarse discretization, are well known in real (and discrete) Josephson ladders, and 
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371 ]. In all our calculations below, is sufficiently large 



have been discussed extensively 
to avoid these spurious steps. 

We discuss first a position- independent coupling between the junction and the cavity. 
An example of the calculated current-voltage characteristics for this case is shown in Fig. 
13] The portion of the IV characteristic for J/Jc < 0.6 corresponds to a soliton which moves 
freely through the long junction. The soliton behaves like a free, massive, but relativistic, 
particle, with limiting velocity c. Its motion is entirely unaffected by coupling to the cavity. 
The free motion of the soliton can be understood in a simple way by the following argument. 
If the cavity electric field is uniform, the coupling between the soliton and the cavity mode 
is independent of the soliton position. Thus, the cavity electric field exerts no force on the 
soliton, which therefore should still travel with constant velocity even if there is a strong 
junction-cavity coupling, consistent with our numerical results. In fact, the IV characteristics 
shown in Fig. E] would be the same if ^ = 0. 

Next, we consider a position- dependent junction-cavity coupling. In this case, the IV 
characteristics are clearly perturbed by the cavity. Figs. |S]|7| shows the soliton branch of the 
full IV curve for a spatially varying coupling of the form used in eqs. fll6p and (fT7j) . By the 
soliton branch, we mean that part of the IV curve produced when the current is increased 
from Jz/Jc = 0.0 to approximately 0.6. We use the same model parameters as for uniform 
coupling, except that we consider several different coupling strengths, namely AE'^g = g', 
with g' = 10~^, 10~^, and 10~^ and the dimensionless cavity frequency Q = 0.75. We also 
take m = 1 in eq. ()19|). There are now clear step- like structures in the IV characteristics for 
all three coupling strengths, at (V) /{NRIc) ~ 0.95, 0.65, and 0.33 respectively, which were 
absent in the case of uniform coupling. These step-like structures correspond in each case 
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FIG. 4: Current-voltage plot for a single long junction in a resonant cavity, with Qj = 10.0, 
A = 0.05, 17 = 0.75, a uniform coupling with qEqA = 1 x 10~^, and Qc = 10. Our initial 
conditions are such that one soliton is present. This curve is the same for any choice of the 
strength parameter qEqS 



to the locking of the soliton to the cavity mode. Specifically, the soliton circulates around 
the long junction at a frequency of one cycle per cycle of the cavity mode. For the weakest 
coupling shown, with g' = 10~'^, the cavity mode is shifted very little from Q = 0.75. The 
voltage step thus occurs approximately at 



NRh Qj 

or (V^) / (NRIc) ~ 0.95 for these parameters, corresponding to one cycle of the soliton around 
the cavity per unperturbed cavity period. For the two stronger couplings shown in Figs. IHl 
and d the soliton is still locked to the cavity, but the coupling is strong enough that the 
cavity frequency is shifted substantially down from its unperturbed value, to about 2f2/3 
and Q/3 respectively. The corresponding time- averaged voltage on the step is approximately 
iy)/{NRh) = {2/3)[Ann/Qj] ~ 0.66, and (l/3)[47rfi/Qj] ~ 0.33 for these two couplings. 

To see the effects of changing the cavity frequency, we have carried out additional calcu- 
lations with the same sinusoidal junction-cavity coupling but different cavity frequencies Q 
and various coupling strengths. The voltage plateaus typically vary approximately linearly 
with Q, as predicted by eq. ()21|). An example of this behavior, for the rather large coupling 
constant g' = 10^^, is shown in Fig. |H1 {V)/{NRIc) is approximately 2/3 the value of Fig. 
as suggested by eq. (j2H). 
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Besides the time-averaged voltages, we have also calculated time-dependent voltage dif- 
ferences at various points across the junction for most of the examples shown in Figs|31IHl 
and others. We have also computed the time-dependent cavity variable Qri^)- In all cases, 
these calculations provide clear evidence of locking between the junction and the cavity. 

Some representative examples of the voltages are shown in Figs. IHl El and ^2 In Figs. 
El^a) and (b), and Figs. [TUT a) and lior b). we show the time-dependent voltages for mini- 
junctions nos. 60, 120, 30 and 90 at 4^ = 0.04, for the parameters and frequency of Fig. |H1 

•J c 

The voltages shown in Figs. Efa) and Efb) have the same characteristic shape and period, 
but are 180° out of phase with one another. The same is true for Figs. El (a) andlTIIlfb). The 
voltages of the second pair differ in wave form, but not in period, from those of the first pair. 
Indeed, we have found that all 120 time- dependent voltages on this step have the same period 
and that the voltages of all mini-junction pairs separated by exactly one-half the junction 
length are identical but 180° out of phase. The corresponding behavior of qr{T) is shown 
in Fig. ^2 Qr{'^) has the same period as that of all the individual mini-junction voltages, 
showing that the long junction is indeed locked to the cavity mode. Behavior similar to that 
of Figs. El and El is also seen in cases with smaller g'. Because of the weaker coupling, the 
amplitude of the cavity parameter on these steps (not shown) is much smaller than in 
Fig. ^2 The period of qr is again the same as that of all the mini-junction voltages, and 
the time-dependent voltages are again identical in pairs, but 180° out of phase, as in Figs. 
HandCni 

We see from Figs. EUHl that the voltage steps are not completely fiat. Despite this 
slight curvature, we have confirmed numerically that the voltages of all the individual mini- 
junctions have the same period as that of the cavity. Furthermore, the time-averaged voltages 
across each mini-junction are all the same. If we examine the voltages at current values off 
the step, we find that the cavity and individual junctions no longer have the same periods. 

The behavior of the time-dependent voltages at different points along the junction can be 
understood from the sinusoidal junction- cavity coupling. For every point on the junction, 
there is a corresponding point separated by Lrc/2 which experiences an equal and opposite 
coupling to the cavity. If the coupling has a spatial dependence sin(27ra;/L^.), there are 
two nodal points along the junction where the coupling is zero, and two points where the 
coupling has maximum amplitude but is 180° out of phase. At these two maximum points, 
the cavity-soliton interaction at any given time is equal in magnitude but opposite in sign. 
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FIG. 5: Current-voltage plot for a single long junction in a cavity, with Qj = 10.0, A = 0.05, 
n = 0.75, N = 120, Qc = 10, and coupling gE^Asm{2Trx/L^), with qEqA = 1. x 10"^. Our initial 
conditions are such that one soliton is present. 




FIG. 6: Same as Fig. El except that gE^A = 1.0 x 10"^. 
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FIG. 7: Same as Fig. El except that gE^A = 1.0 x 10"^^ 
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FIG. 8: Same as Fig. [71 except that = 0.5. 
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FIG. 9: Time-dependent voltage V{t) for small junctions no. (a) 60 and (b) 120, for 4000 < r < 
4050, plotted versus time r. Current is such that time-averaged voltage is on the lowest step-like 
structure of Fig. 8. 



Because of this feature, the time-dependent voltages at these two points should have the 
same wave form but should be 180° out of phase, as is seen in Figs. M (a) and (b). 

This picture also explains why the time-dependent voltages along the junction are at 
all points equal in pairs but 180° out of phase. The different pairs have different voltage 
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FIG. 10: Same as Figs. El (a) and (b) except that we plot voltages for junctions no. (a) 30 and (b) 
90. 
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FIG. 11: Cavity variable Qrir) for the long junction of Fig|Hl plotted versus r for the voltage on 
the lowest step-like structure for 4000 < r < 4050. 
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FIG. 12: Current-voltage plot for a single long junction in a cavity, with same parameters as in Fig. 
|71 except that Cl = 1.5 and the coupling is ^E'q A sin(2/cx) with k = 2tt/Lx and gE^A. = 1. x 10^^. 
Our initial conditions are again such that one soliton is present. 

wave-forms because the coupling amplitude between each mini-junction and the cavity varies 
with the spatial dependence sin(fcx). Despite the different wave-forms, we have verified that 
the time- averaged voltage difference is the same at each point along the junction. This is 
consistent with the picture that this voltage is produced by a soliton which passes each point 
along the junction with the same average frequency. 

Further evidence of a strong soliton-cavity coupling can be seen by comparing Figs. 
[7|and^l Fig. [7[ as noted above, shows the soliton branch of the IV curve for Q = 0.75, 
qEqA = 10~^, and a coupUng varying spatially as sm{27ix/Lx). In the second, we assume the 



18 



same coupling strength of qEqA = 10~^ but a frequency Q = 1.5 and a spatial dependence 
of sin(47ra;/L^.). In both cases, there are steps at approximately the same value of {V), 
namely {V)/{NRIc) ~ 0.33 ~ vr/lO. These figures show that simultaneously doubling the 
cavity frequency and halving the coupling periodicity leaves the step height unchanged. We 
interpret this behavior as showing that, in both cases, the soliton is locked to the cavity 
mode so that it moves by a distance equal to the wavelength of the mode per mode cycle. 

We have also computed qr{T) for J^/Jc = 0.05 and the parameters of Fig. 1121 fnot shown). 
We find that this amplitude is much smaller than that shown in Fig. 11 for the step of Fig. [7| 
A simple qualitative explanation for this behavior is the following. The soliton has the same 
spatial width in both cases, but the junction-cavity coupling varies spatially more rapidly 
in ^1 than in Fig. [7| Since the junction-cavity coupling thus varies substantially over the 
width of the soliton in Fig. ^[ it has a smaller effect than in Fig. [7| Thus, we expect a 
much weaker step-like feature in the IV characteristic of Fig. than in Fig. [71 as is indeed 
observed numerically. 

IV. DISCUSSION 

A. Qualitative Argument for Soliton Steps 

The locking of the moving soliton to the mode of the resonant cavity can be accounted 
for by a simple analytical argument. The argument starts from the equations of motion ()16|) 
and (fT7|) . We assume solutions of the form 

t) = 00 + k{x - vt) + 01 sin[A;(a; - vt)] (22) 
qr{t) = Re(goe'"*) (23) 

This solution is suggested by the fact that, if g' = 0, the solution 4>{x,t) is rigorously of the 
form (f){x,t) = (f){x — vt), where v is the soliton velocity. We now substitute these assumed 
solutions into eqs. (fTBj) and (fT7|) . carry out the indicated derivatives, and use the standard 
expansion 

oo 

sin {(pQ + k{x — vt) + (j)ism[k{x — vt)]} = ^ Jn{(pi) sm[(j)Q + k{x — vt)+nk{x — vt)], (24) 
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where Jn(0i) is the n*'* order Bessel function. In eq. (fT7|). we use the expression (pu = 

—uj'^ipi sm[k{x — ft)], where uj = vk, and carry out the integral to obtain 

Qr + -p^Qr + ^"^Qr = COs(u;t). (25) 

Since this equation is hnear, the driven solution for is simply 

Solution (PU)) can now be substituted, along with the Bessel function expansion, eq. (j21I), 
back into the other equation of motion, eq. (fT!)|l . Next, we assume that the amplitude 0i is 
small, and expand the Bessel functions in powers of 0i. Finally, we retain only the lowest 
Fourier components, namely those involving a constant term, and the functions sm[k{x — vt)] 
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and cos[k{x — vt)]. A similar approach has been used previously in Refs. 
treat intrinsic Josephson junctions to optical phonons in high-temperature superconductors. 
Setting the coefficients of each of these terms separately equal to zero in the expanded 
version of eq. (fT^ . we finally obtain the following three equations for 0o, 0i, and (J/Jc): 

1 K' 

-J sin 00 + '77.''^'' = 0; (28) 



cjpcj , 1 , K"iVLuj/Q 



-k\l- ^)0i + 4 cos 00 + ^ ' ^ = 0, (29) 



where u = vk, k = 2n/L^, K" = {h^c'^g'yuj^Lj[167rd{e*yM] and 

D = {Q^- u^f + nW/Ql (30) 



is a resonant denominator |41|. 

Eqs. fl28|) and ()29|) are readily solved for 0o aiid 0i in terms of the soliton velocity v, and 
the result substituted into eq. (j2ZI)- The resulting eq. (P7j) expresses the current J in terms 
of the soliton velocity v, or equivalently, the time-averaged voltage across the long junction, 
{V)t = 2nhv / {2eLx) . As is suggested by the resonant form of eqs. and (jSHI), this IV 
characteristic has a peak when u ^ Q, corresponding to the voltage plateaus observed in 
our numerical results. 

This same approach also shows, through eq. ()26p. that there is a peak in the amplitude 
Qr of the cavity mode when the same condition (a; = Q) is satisfied. Once again, this peak 
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in the amplitude is observed in our numerical simulations. For example, Fig. 11 shows the 
large amplitude of qr at a voltage satisfying the resonance condition; the amplitude of at 
other voltages is much smaller. 



B. Comparison to Predictions of Other Models 

It is useful to compare our model and results for the coupled soliton-cavity system to 



those of other workers. For the case o: 
to those obtained in Refs. UM and 



a^uniform coupling, our results are formally analogous 
These workers consider the response of intrinsic 
Josephson junctions in a high-Tc superconductor coupled to an optical phononmode within 
the junction. For spatially varying coupling, our numerical results for the locking of a 
moving fiuxon to a cavity mode somewhat resemble those of Refs. |2S] and [2^, obtained 
for a coupled fluxon/optical phonon system in an intrinsic Josephson junction, though the 



equations describing the two systems are not identical 



However, there is a significant 



difference in the physics. The anomalies found in Refs. Q and Q: 

in these systems, the 

quantities which play the role of cavity modes are the optical phonons, which are intrinsic 
to the junctions themselves. By contrast, our cavity modes are assumed to arise from some 
cavity extrinsic to the junctions. 



C. Possible Realizations of Sinusoidal Coupling 

Finally, we briefly discuss the type of electromagnetic mode which could produce the 
sinusoidal coupling we use. The electric field of this mode has a non-zero curl which varies 
sinusoidally with position, i. e., it has a sinusoidally varying magnetic field. This type of 
mode should be readily achievable in a real cavity. The voltage steps should be achievable so 
long as the phase velocity in the cavity is smaller than the limiting soliton velocity c. Another 
way to produce this type of coupling would be actually to prepare a long Josephson junction 
in the shape of a ring, and then to use a cavity mode with a spatially uniform magnetic 
field. In this case, the junction occupies the ribbon-like region between two circular rings of 
superconductor; if the planes of the rings are parallel to the z axis, then the magnetic field 
of the cavity mode should be uniform and parallel to one of the ring diameters. This will 
produce a flux through the long junction which varies sinusoidally around the ring, as in 
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our model. Measurements using a ring geometry, and a static uniform magnetic field, have 



who also show a schematic 



recently been reported, in another context, by Wallraff et al. 
picture of this geometry. 

V. SUMMARY 



In this paper, we have described a model for a long underdamped Josephson junction 
interacting with a single-mode electromagnetic cavity. In our model, we have assumed a 
capacitive interaction between the junction currents and cavity mode, but we consider both 
a uniform coupling and one which varies spatially along the junction length. If no soliton (i. 



e. no fluxon) is present, the junction behaves very much like a small Josephson iunction|ll|. 
In particular, there are SIRS's just as in a small junction, which occur at the voltages 
expected for a small junction. 

If a soliton is present, and the junction-cavity coupling is position- independent, then we 
find that there are no SIRS's. The absence of SIRS's in this case is easily understood: since 
the soliton-cavity interaction energy is independent of the soliton's position, the cavity exerts 
no force on the soliton. If, however, the coupling varies sinusoidally with position, we find 
step-like structure in the IV characteristics, which arise from the junction- cavity coupling. 
These structures arise from the locking of the soliton to the cavity mode, so that the soliton 
travels a distance equal to one wavelength of the coupling interaction during one cycle of 
the cavity mode, or equivalently, the average soliton velocity equals the phase velocity u/k 
of the cavity mode. We find clear evidence of the locking between the cavity mode and the 
junction through the time-dependent voltages across the junction. Everywhere along the 
junction, these time-dependent voltages have the same period of oscillation as the cavity 
mode. Moreover, the voltage differences of points separated by half a wavelength have 
exactly the same wave form, but are 180° out of phase, showing that the soliton travels one 
wavelength per cycle of the cavity mode. 

We have also presented a simple qualitative argument which explains both the positions 
of the self-induced voltage steps and the occurrence of a peak in the radiated energy on the 
steps. This argument agrees well with our numerical results in the limit of weak junction- 
cavity coupling. 

If the junction-cavity coupling is strong, this simple argument does not give the position 
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of the anomalies in the IV characteristics. Instead, the voltage of the step-like structure 
is shifted substantially down, and the frequency of the cavity mode is strongly red-shifted. 
Nonetheless, the cavity is still locked to the soliton motion. 

The voltage on the steps is not absolutely constant, but varies slightly with current. 
When the current lies on a voltage step, the cavity is strongly excited, with a large time- 
averaged squared amplitude g^. If the current does not lie on a step, cavity and Josephson 
junction are not locked, and is much smaller (at least by an order of magnitude in all our 
numerical runs). This behavior is once again in agreement with our simple analytical model 
of the previous section. 

The present model could be modified to apply to a stack of long junctions coupled to 
a single-mode cavity. In the absence of a cavity, it has been known for some time that 
very anisotropic high-Tc cuprate superconductors, such as BiSr2Ca2Cu208+x, behave like 



a stack of underdamped Josephson junctions 
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Coupling such a stack to a cavity is of 



great interest because it may provide a means for phase-locking these junctions, and hence 
providing a coherent source of sub-THz radiation. The dynamics of a stack of long junctions 
without a cavity has previously been modeled by Sakai et al 40]. A stack of long junctions 
coupled to internal phonon modes within the junctions has been modeled by Preis et a^^. 
The present work suggests a means of modeling the coupling of a stack of junctions to the 
same electromagnetic cavity. 
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VII. APPENDIX: ALTERNATE DERIVATION OF LAGRANGIAN FOR LONG 
JUNCTION COUPLED TO A CAVITY. 

In this Appendix, we present an alternate derivation of the cavity-junction Lagrangian 
obtained in Section II. The derivation here is more general, in that we consider coupling to 
both the electric and magnetic fields of the cavity. 
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We write the total Lagrangian as 

L = - Ui. (31) 

The kinetic energy, K, is written as 

A-,=/<i'x!(^. (32) 

where e(x) is the (possibly position-dependent) dielectric function, and E(x) is the electric 
field. We assume that the electric field is the sum of two parts: that due to the junction, 
which we denote Ej„„c, and the part due to cavity mode, which we write as Ecav The 
junction field takes the form 

Then takes the form 



Here the volume integral is taken over the cavity, and we are assuming a geometry such that 
the junction is contained within the cavity. 

The potential energy Ui is the sum of two terms. The first is the Josephson energy, 

C/i,j = -Ly j cos 0. (34) 

The other part is the energy of the magnetic field. This may be written 

Ul,mag = / ^''x^^. (35) 

Once again, the magnetic field may be written 

where ^junc a-nd ^cav a-re the fields due to the junction and cavity. As in the text, we again 
assume that J„ = J^, so that 

^^"'^^ = ^^^^ 
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Then the magnetic field energy takes the form 



2^2 



1 $0 



To make further progress, we introduce operators describing the fields ^cav and 
Both may be expressed in terms of the operator for the cavity vector potential: 



Aea.(x,t) = ( ^ J [a(t) +at(t)]E(x). 



(38) 



Here a and are the annihilation and creation operators for the cavity mode, which satisfy 

the usual Bose commutation relations 

[a, a+] = 1. (39) 

E(x) is proportional to the position-dependent electric field of the cavity mode; its normal- 
ization is given below. Q is the frequency of the cavity mode. In terms of this operator, the 
electric field operator Ecai;(x, i) is given by 

EU^,t) = = -i{hny/'{a - at)E(x) (40) 



and the magnetic field operator Bca^(x, i) is 



Bea.(x, i) = V X A,„„ = l^—j [a + at]V x E(x). (41) 
The operator describing the total energy in the cavity is 



cav 



d X , (42) 



We calculate the ensemble average of this operator, using the results {a^a) = n; (aa^) = n+l; 
(aa) — (a^a^) = 0, with the result 



|E|2 + — |V X Ep 



d^x. (43) 



In order for this energy to equal hil{n + |), we require that the function E(x) be normalized 
so that 

J |E|2 + ^|VxE|2 d^x^2. (44) 
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Also, we assume that the cavity energy is equally distributed between the electric and 
magnetic fields. This imphes that 

j \E\'^d^x = ^ y |V X El^d^x = 1. (45) 

Having obtained the operator forms for the cavity electric and magnetic fields, we are 
now in a position to derive the electric and magnetic parts of the junction-cavity coupling. 
After some algebra, the electric field part of the coupling Lagrangian may be written 

Lcoup,E = J dx^(t>t[-i{hQy^^a - a'')]E,{x). (46) 

Similarly, the magnetic field part of the coupling Lagrangian may be written 

Lcoup,B = ^Ly^ihcyny/^ J dx{a + at)0,(V x E(x)),. (47) 

To make contact with the notation in the main part of this paper, we introduce position 
and momentum operators for the cavity mode, by 

Vr = l^-^j ^(at - a) (48) 

and 

''=(2^0) + 

These operators have the standard canonical commutation relations 

[Pr, Qr] = -i^- (50) 
In terms of these operators, we may write the Hamiltonian for the cavity alone as 

Hcav J (Ecai) ■ ^cav ~l~ ^cav ' ^cav^ 2,M ~^ 2 ' (^'^) 

where we define the "spring constant" 

K = Mfl^. (52) 

Prom Hcay we may infer the corresponding cavity Lagrangian using — Pr/M, with the 
result 

Lcav = ^ - l^ql (53) 
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It is readily verified that classical Hamiltonian equations of motion resulting from the Hamil- 
tonian (j^Tj) are equivalent to the classical Lagrange equation of motion obtained from the 
Lagrangian both lead to 

Qr - n'^Qr = 0. (54) 

Finally, we can express the coupling Lagrangians in terms of the variables introduced 
above. First Lcoup,E takes the form 

Lcoup,E ~ J ^•'^^coup,e{,x)i (55) 

where 

C-coup,e{x) = ~g{x)qr(f>t, (56) 

and 

n(T^ = —f\ / 

47r2e 



(7(x) = -eW^|-E(x)-z. (57) 



Similarly, Lcoup,B takes the form 

Lcoup,B Ly / dxCcoup,B{x), (58) 

where 

C-coup,b{x) = -gB{x)qr(l)x, (59) 

and 

gBix) = ^^%VxE-y. (60) 

There is one additional term in the Lagrangian which also represents a coupling between 
the cavity and the junction. This is the term 

f d 

Ld = Ly J dx — {e - l)Ecav ■ Ecav (61) 

But this term does not couple the cavity variables to those of the junction. Instead, its only 
effect will be to produce a slight shift in the cavity resonance frequency. We have therefore 
not included this term in our calculations described in the body of the paper. 
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